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ABSTRACT 


The response of a composite cylinder (metallic cylinder coated 
with a rubber material) subjected to an underwater explosion was 
analyzed numerically. Qualitative differences between coated and 
uncoated cylinders were investigated. The dynamic response of the 
coated cylinder was found to be adversely affected when impacted by 
an underwater shock wave under certain conditions of geometry and 
material properties of the coating. When adversely affected, 
significant deviations in axial and hoop stress and strain values 
were observed. The coated cylinder exhibited larger effective 
plastic strain and higher residual internal energy in the metallic 
material. Rubber coatings appeared to inhibit energy dissipation 
from the metallic material to the surrounding water medium. A 
parametric study of various coatings was performed on both aluminum 
and steel cylinders. The adverse effects of the coating decreased 
when the shear modulus of the rubber was increased or when the 
rubber thickness was increased, indicating the existence of 
threshold values for these parameters. The results of this study 
indicate that the stiffness of the coating is a critical factor 
involving these threshold values. 
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I. INTRODUCTION 


Not since World War II has mankind witnessed such dramatic 
and profound changes as those that have occurred in recent 
years. Sweeping political and economic reform in the former 
Soviet Union has not only slowed their military production and 
affected their deployment of forces, but it has been the 
impetus that has led the United States to re-evaluate its own 
strategic planning and policies. Accordingly, the United 
States Navi’' has had to re-formulate its “Maritime Strategy". 

Although the Soviet threat has diminished, the United 
States Navy must remain capable of quickly responding to 
conflicts around the world, as demonstrated in Operation 
Desert Shield and Operation Desert Storm. To continue to meet 
all its commitments abroad, amidst a shrinking defense budget, 
the Navy has placed a premium on development and production of 
ships and submarines that are technologically superior, cost- 
effective, easily maintained and perhaps most importantly, 
survivable. In the 21st century, "shallow water" naval 
engagements are most likely. In this environment, U.S. 
warships and submarines will be most vulnerable to mines and 
depth charges. Although a future adversary may not possess a 
large arsenal of military firepower, it will likely employ 
sophisticated weapons, possibly nuclear weapons, capable of 
inflicting considerable damage without rendering a direct hit. 


1 






To counter this threat, the United States Navy has 
intensified its efforts in research and study of underwater 
explosions on ship and submarine hulls. No matter how well a 
ship or submarine might be armored, the effects of a direct 
hit from a mine, torpedo or depth charge are potentially 
catastrophic. If however, the explosion were to occur away 
from the hull, then the ability of the hull to withstand the 
effects of the shock wave caused by the blast becomes crucial. 
In an attempt to better understand the effects of blast damage 
on shock hardening, the Navy has performed underwater 
explosion tests and analyzed the problem numerically using 
computer simulations. Knowledge gained through such study 
should prove invaluable to ship and submarine designers of the 
future. 

The dynamic response of a structure to an underwater 
explosion is a complex problem due to the interaction between 
the fluid and structure. Over the past several years, the 
Naval Postgraduate School has conducted extensive research in 
the field of underwater shock. This study has involved both 
physical testing and numerical modeling of uncoated metal 
cylinders of various configurations such as unstiffened, 
stiffened, single shell layer and double shell layer models. 
Ship shock qualifications and even small scale testing can be 
cost prohibitive and time consuming, therefore much of the 
work has focused on developing computer models using the 
numerical analysis technique. Results from earlier studies 
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[Refs. 1-3] indicate response predictions compare favorably 
with both analytical solutions and experimental data. 

In a recent study conducted by the Navy [Ref. 4], when 
subjected to an underwater shock wave, damage to a steel test 
panel increased when it was covered with a low density layer 
of glass microspheres suspended in water (slurry). This result 
was unexpected and the reason for this adverse response is 
still under investigation. It is suspected that compliant 
coatings may affect the propagation of a shock wave and the 
dynamic response of a structure. 

The objective of this study was to examine the response of 
a metal cylinder coated with a rubber material, when subjected 
to an undei'water explosion, utilizing the numerical analysis 
technique. Rubber coated aluminum and steel cylinders were 
analyzed and a parametric study of various coatings was 
performed to gain a better understanding of the damage 
mechanism. Preliminary* study seems to indicate that the 
dissipation of energ:^" from, the structure to the surrounding 
water medium is critical to its dynamic response. 
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II. NUMERICAL ANALYSIS 


A. NUMERICAL CODES 

To Study the effects of an underwater shock wave on a 
cylinder, a coupled finite element and boundary element code 
was utilized. The two codes operate in tandem so that the 
fluid-structure interaction can be calculated. This code was 
successfully implemented at the Naval Postgraduate School in 
1991 and has since been the primary tool used for studying 
underwater shock phenomena. Structural calculations were 
handled by the finite element code while the effects of the 
propagation of the explosive pressure wave through the water 
medium were calculated by the boundary element code. 

1. Finite Element Code 

The finite element code used in this analysis was 
VEC/'DYNA3D [Ref. 5]. Originally implemented in 1976, it has 
been widely used in the industrial sector and at many 
universities. The code is efficient and extremely flexible, 
offering a wide variety of material models from which to 
choose. 

Input to VEC/DyNA3D was provided by LS-INGRID [Ref. 
6], a pre-processor and three dimensional mesh generator. 
Model configuration, element and material type, boundary 
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constraints and loading conditions were specified in the pre¬ 
processor . 

Output from VEC/DYNA3D was provided by LS-TAURUS [Ref. 
7], an interactive post-processor that was used to display 
fringe plots, element, node and material time history plots, 
and generate data files. 

2. Boundary Eleaent Code 

The boundary element method code used in this analysis 
was USA (Underwater Shock Analysis) (Ref. 8]. The USA code 
computed the transient response of the submerged cylinder to 
the shock wave. Calculations of the fluid-structure 
interaction are based on the Doubly Asymptotic Approximation 
(DAA) theory developed by T.L. Geers [Ref. 9]. The boundary 
element method precludes the requirement to discretize the 
water medium since equivalent forces and masses are 
transferred to the cylinder from the water and are applied at 
the nodes of a two dimensional mesh which is superimposed on 
the cylinder surface. As a result, computational efficiency is 
improved. 

3. Structure and Wave Equations 

The differential equation for the dynamic response of 
a structure can be expressed as follows: 

[M,] bei + [c,] ui + [if,] ui = if) (1) 

w’here [/fj , [c,] and [if,] are the structural mass, damping and 


5 





stiffness matrices respectively, and , biJ and W are the 
nodal acceleration, velocity and displacement vectors 
respectively. If} is the excitation force vector and can be 
expressed as a function of the incident and scattered 
pressures of the impinging shock wave and any concentrated 
loads applied to the structure. Their functional relationship 
is as follows: 


{f} = - [G] [Af] {Pi* PJ ^ {fj (2) 

where [G] is the fluid/structure transformation 
matrix, [Af] is the diagonal area matrix associated with the 

fluid elements, iP^+P^ is the incident and scattered pressure 

wave vector and is the force vector applied to the dry 

structure. The scattered pressure represents the only unknown 
quantity, but it can be determined by solving the first order 
wave equation given below, 

[Mt] {PJ + pc[Af] {Pj = pc[Wf3 iiij (3) 

where , p , c and (U^ are the symmetric fluid mass matrix, 

fluid density, sound speed and particle velocity vector 
respectively. Although the general solution is complicated, 
early and late time solutions to equation (3) can be found by 
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approximating {p^> >{Pg} and {Pg}< <{Pg} respectively and are 

referred to as the Doubly Asymptotic Approximations. The early 
time approximation is also referred to as the high frequency 
approximation or plane wave approximation. The solution is 
calculated at early times, at locations close to the structure 
and assumes the shape of the shock wave to be planar. The 
change in the scattered wave pressure is almost instantaneous, 
so becomes the dominant term. Simplification and 

integration of equation (3) leads to the following solution: 

{pj = pciUj,) (4) 

The late time approximation is also referred to as the lew 
frequency approximation or virtual mass approximation. The 
solution is calculated when the scattered pressure wave has 
travelled a significant distance from the structural surface 
so that the pressure change becomes negligible and 
the IP^ term can be ignored. Simplification of equation (3) 

leads to the following form: 


[A,] {Pj = [Mf] iii; (5) 

Determination of intermediate time solutions is accomplished 
by using the DAA relations which bridge the upper and lower 
limit approximations. A second DAA theory, called DAA2 was 
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later developed to account for surface curvature effects. DAA2 
theory is contained in [Ref. 10]. 

B. NUMERICAL MODEL 

1. Material Models 

The behavior of the metal material was based on the 
Kinematic/Isotropic Elastic-Plastic Model [Refs. 5,6]. In 
addition to the density, the model required at least two other 
material properties be specified such as the bulk modulus. 
Young's modulus, shear modulus or Poisson's ratio. Optional 
inputs included yield stress, hardening modulus, or hardening 
parameter. The hardening parameter reflects the type of 
hardening (Isotropic or Kinematic) the material will 
experience. When steel was analyzed, strain-rate hardening 
effects were accounted for and were based on the Cowpers- 
Symonds Model, shown in equation (6) 

( 6 ) 

where Oy is the yield stress and c and p are the strain-rate 

parameters which were obtained from a previous study (Ref. 3] . 

The behavior of the rubber material was based on the 
Compressible Mooney-Rivlin Model [Refs. 5,6]. Mooney pioneered 
a new approach to analyze the deformation of a soft material 
such as rubber [Ref. 11]. He stated that classical elastic 
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theory could not be applied to a highly elastic (hyper-elastic 
or superelastic) material but that deformation could be 
accurately represented in terms of its general strain energy 
density. 

Mooney postulated that in addition to being homogeneous 
and free from hysteresis, a hyper-elastic material possesses 
the following properties: 

It is isotropic in the undeformed state and remains 
isotropic in planes that are at right angles to a stretch 
or squeeze. 

Deformations are isometric, i.e. occur without change in 
volume. 

- In simple shear, shearing stress in any isotropic plane is 
proportional to the shear force. 

Mooney expressed the general strain energy* density function of 
a material, Was follows: 








-3) 


(7) 


where Xj is the principal stretch, which is equivalent 

to (1*€|) , where is strain. Applying conditions for a 

hyper-elastic material and substituting expressions 
for Cj and C, as shown below 




(G+ / i 


4 


. _ iG-H) 
a A 


( 0 )*< 9 ) 
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leads to the following form: 



( 10 ) 


where G and H represent the modulus of rigidity and the 
modulus characterizing asymmetry of reciprocal deformation 
respectively. The variable H is a measure of the ability of a 
material to store energy when compressed as opposed to when 
stretched. To express the asymmetry of reciprocal deformation 
in a more useful form, Mooney defined a new parameter a , 
called the coefficient of symmetry, as follows: 


a- 


H 

G 


( 11 ) 


Experimental values for G ar i a w'ere derived for two types of 
rubber, tread stock and gum stock. Mooney found that the 
experimental data relating force and elongation, agreed 
closely with analyt:.':al results for soft rubber undergoing up 
to 400"o elongation and 50% compression. Mooney believed that 
just as the infinitesimal deformation of a hard cr moderately 
elastic material can be characterizad by two material 
properties, shear modulus and Poisson's ratio, deformation of 
rubber can also be characterized by two properties, shear 
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modulus and coefficient of asymmetry. Mooney's analysis and 
development form the basis upon which the numerical model for 
the rubber material was generated. In the finite element 
model, the strain energy density function was defined as 
follows. 

W= A{I-3) + B(II-3) + C(JIJ-2-l) +15(1X1-1)2 (12) 


where 

C=0.5A+B , 7^ I , (13).(14) 

2 {l--2v) 

2 (A+B) - G (15) 

V is Poisson's ratio and are invariants of the right 

Cauchy-Green tensor. The model required values for the 
constants A and B . By relating the first two terms of the 
strain energy' density equation, equation (12), with equation 
(7), as expressed by Mooney, equivalent expressions for the 
two constants were determined and are shown below. 

A=^(l+a) . B=-?(l-a) (16). (17) 

4 4 

Since value? for G and o were given for tread stock and gum 
stock in Mooney's development, these types of rubber were used 
in the analysis. 
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2. Simulated Tei£t Conditions 

All cases involved identical test conditions. The 
problem simulated detonation of 60 lbs. of HBX-1 spherically- 
shaped explosive at a horizontal standoff distance of 25 ft. 
from the cylinder resulting in a side-on attack geometry. The 
uncoated cylinder models (metal only) were 12 in. in diameter, 
42 in. long, 1/4 in. thick (shell thickness) with a 1 in. 
thick endplate. The composite or coated cylinder models had 
the same dimensions for the metal material as the uncoated 
cylinders but in addition were coated with 1/4 in. of rubber 
material on the outer surface of the cylinder. The shell and 
coating thicknesses were nominal values and were changed for 
parametric study. The cylinder geometry/test profile is shown 
in Figure 1. 



rigur« 1: Cylinder geometry and test profile. 
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Actual pressure-time history data from a previous physical 
test was utilized in the analysis.' The pressure-time history 
plot is shown in Figure 2. 


FREE-FIELD PRESSURE-TIME HISTORY 



Figure 2: Free-Field Pressure-Time History Plot. 

3. Finite Element Mesh 

The problem geometry involved two planes of symmetry 
(x-z and y-z planes) which allowed the cylinder to be modeled 
as a quarter cylinder. This significantly reduced the number 
of computations in each test. Appropriate boundary constraints 
were applied along the symmetric planes. The finite element 
mesh for the quarter model is shown in Figure 3. The metal 
material consisted of 186 thick shell elements while the 
rubber material consisted of 198 solid (brick) elements. 
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Figure 3; Finite element mesh (quarter model). 




Only solid elements could be used for the rubber in this model 
configuration. A total of 9 locations were analyzed in each 
case, 3 at the center section of the cylinder, 3 near the 
endplate and 3 midway between the center section and the 
endplate. At each section the 3 locations were approximately 
equally spaced around the surface of the cylinder. A1 
represents the element closest to the charge. The metal shell 
material in Figure 3 shows the element orientation. Unless 
otherwise specified, all plots reference these element 
locations. Although the response of the rubber was examined in 
different cases, the analysis of the metal response was the 
primary objective, therefore most of the results presented are 
for the metal material. 

4. Element Compatablllty 

Since the metal and rubber were modeled using a 
different type of element, it was necessary to verify the 
elements were compatable and would give reliable results. This 
was done by performing a test simulation for three different 
shell layer configurations, a single thick shell model, a 
double thick shell (two adjacent single thick shells) model, 
and a thick shell and brick model, then comparing the 
responses at various locations. The total shell thickness was 
the same in all configurations. Results are given in Chapter 
III. 
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5. Material Properties 


Nominal material properties were obtained from a 
materials handbook. Aluminum, steel and rubber material 
properties are shown in Tables 1-3 respectively. 


TABLE 1: ALUMINUM (6061-T6) PROPERTIES 


Parameter 

Symbol 

Value 

Density 

P 

5.412 slugs/ft^ 

Poisson's ratio 

V 

0.33 

Young's modulus 

E 

1.08x10’ psi 

Yield stress 


4.0x10'* psi 

Speed of sound 

1 ^ 

16,400 ft/sec 


TABLE 2: STEEL (ASTM A106 grade C) PROPERTIES 


Parameter 

Symbol 

Value 

Density 

P 

15.218 slugs/ft^ 

Poisson's ratio 

V 

0.30 

Young's modulus 

E 

2.9x10’ psi 

Hardening modulus 

^tan 

1.114x10^ psi 

Yield stress 


4.7x10'* psi 

Strain-rate parameter 

c 

8.797x10^ sec'^ 

Strain-rate parameter 

p 

5.65 

Speed of sound 

a 

16,900 ft/sec 
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TABLE 3: RUBBER PROPERTIES 


Parameter 

Symbol 

Value 

Density 

P 

1.908slugs/ft^ 

Poisson's ratio 

V 

0.4991 

Speed of sound 

a 

100 ft/sec 

Shear modulus (Tread stock) 

G 

95.8 psi 

Coef. of asymmetry (Tread stock) 

a 

0.223 

Shear modulus (Gum stock) 

G 

48.6 psi 

Coef. of asymmetry (Gum stock) 

a 

0.448 
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III. RESULTS AND ANALYSIS 


As discussed in Chapter 2, data was collected at nine 
different locations on the cylinder, corresponding to the 
element positions as shown in Figure 3. Results are reported 
for both axial and hoop directions. The axial direction is 
parallel to the centerline axis of the cylinder while the hoop 
direction is tangent to the circumference of the cylindrical 
shell for a given location. Some data has not been reported 
since the results are often redundant and can be seen in other 
plots. For this reason, strain plots are shown in lieu of 
stress plots since the same features and trends are exhibited 
in both. 

A. ELEMENT COMPATABILITY 

A preliminary study was undertaken to verify element 
compatability since the metallic material and the coating were 
modeled using different types of elements. Single thick shell, 
double thick shell and thick shell/brick element 
configurations were compared at all "A" and “B" locations and 
are shown in Figures 4-15. Axial and hoop strain values for 
the single thick shell and double thick shell cases were 
essentially identical. The thick shell/brick model deviated 
slightly from the thick shell models in the magnitude of the 
response. However, their frequency of oscillation appeared to 
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AXIAL STRAIN FOR ELEMENT A1 



TIME(sec) 

rigure 4: Axial strain at position A1 for different shell 
thickness element configurations. 


HOOP STRAIN FOR ELEMENT A1 



Figure 5: Hoop strain at position A1 for different shell 
thickness element configurations. 


19 
























































MICROSTRAIN 


AXIAL STRAIN FOR ELEMENT A2 



Figure 6; Axial strain at position A2 for different shell 
thickness element configurations. 


HOOP STRAIN FOR ELEMENT A2 



TIME(s«c) 

Figure 7: Hoop strain at position A2 for different shell 
thickness element configurations. 
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AXIAL STRAIN FOR ELEMENT A3 



8; Axial strain at position A3 for different shell 
thickness element configurations. 


HOOP STRAIN FOR ELEMENT A3 



TIME(»m) 

Figure 9: Hoop strain at position A3 for different shell 
thickness element configurations. 
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Figure 10; Axial strain at position B1 for different shell 
thickness element configurations. 


HOOP STRAIN FOR ELEMENT B1 



TtME(Me) 

Figure lit Hoop strain at position B1 for different shell 
thickness slernent configurations. 
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rigurtt 13: Hoop strain at position B2 for different shell 
thickness element configurations. 
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be approximately the same. At all locations, the thick 
shell/brick model appeared to approach the steady state more 
rapidly. The results indicate there was no interfacing problem 
or unusual behavior caused by the two different types of 
elements and demonstrated that the thick shell/brick model 
could be expected to behave consistently and provide reliable 
response data. The thick shell and brick (solid) elements were 
therefore considered to be fully compatable and were used in 
the subsequent study. 

B. EFFECT OF StTRFACE COATINGS ON CYLINDERS 

The purpose of this study was to determine the effects of 
surface coatings on the dynamic response of cylinders 
subjected to underwater shock. A total of six cases were 
analysed. In the first case, an uncoated aluminum cylinder was 
compared to two composite (coated) cy.linders, one coated with 
tread stock rubber (referred to as Composite Model 1) and the 
other coated with gum stock rubber (referred to as Composite 
Model 2). The aluminum shell material and the rubber coating 
were both 0.250 inches thick. Both composite cylinders 
exhibited higher values of stress and strain than the uncoated 
cylinder. In early time, (less than 1 msec.) response was 
approximately the same. (Passage of the shock wave occured at 
about 0.2 msec.). At most locations, the deviations became 
significant after 1 msec.. Axial and hoop strain were plotted 
for all elements of interest and are shown in Figures 16-33. 
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AXIAL STRAIN FOR aEMENT A1 



Piffure 16: Axial strain at position A1 for uncoated and coated 
alurr.inuxa cylinders. 


HOOP STRAIN FOR ELEMENT A1 



Pigfur« 17: Koop strain at position A1 for uncoated and coated 
aluminum cylinders. 
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Figure 18; Axial strain at position A2 for uncoated and coated 
aluminum cylinders. 


HOOP STRAIN FOR ELEMENT A2 



Figure 19: Hoop strain at position A2 for uncoated and coated 
aluminun^ cylinders. 
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AXIAL STRAIN FOR ELEMENT A3 



TIME(sec) 

Figure 20; Axial strain at position A3 for uncoated and coated 
aluminum cylinders. 


HOOP STRAIN FOR ELEMENT A3 
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Figure 21: Hoop strain at position A3 for uncoated and coated 
aluminum cylinders. 
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AXIAL STRAIN FOR ELEMENT B2 



Figure 24: 


Axial strain at position B2 for uncoated and coated 
aluminum cylinders. 


HOOP STRAIN FOR ELEMENT B2 



TIME(ssc) 


Figure 25: Hoop strain at position B2 for uncoated and coated 
aluminum cylinders. 
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AXIAL STRAIN FOR ELEMENT B3 



figure 26: 


Axial strain at position B3 for uncoated and coated 
aluminum cylinders. 


HOOP STRAIN FOR ELEMENT 83 



TIME{»c) 

Figure 27: Hoop strain at position B3 for uncoated and coated 
aluminum cylinders. 
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Figure 28; Axial strain at position Cl for uncoated and 
coated aluminum cylinders. 
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Figure 29: Hoop strain at position Cl for uncoated and coated 
aluminum cylinders. 
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Figure 30; Axial strain at position C2 for uncoated and coated 
aluminum cylinders. 


HOOP STRAIN FOR ELEMENT C2 


ALUMINUM MODEL 
COMPOSITE MODEL 1 
COMPOSITE MODEL 2 




TIME(tftC) 

Figure 31; Hoop strain at position C2 for uncoated and coated 
aluminum cylinders. 
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Neither composite cylinder showed consistently higher 
magnitudes than the other, as this varied depending on the 
location analyzed. Both composite models exhibited higher 
oscillation than the uncoated model. The results from this 
study seemed to indicate that the rubber coating tended to 
concentrate energy into the metal cylinder rather than allow 
for a rapid release of energy into the water mediumi as in the 
uncoated case. The effects of the rubber coating are clearly 
visible in the deformation the cylinders. Figures 34 and 35 
show deformation of the uncoated cylinder and composite 
cylinders 1 and 2 at 2.86 msec, and 5.86 msec, respectively. 
The unusually high values of strain attained at locations Cl 
and C2 may be attributed to the inertial effects of the thick 
endplate. Figure 36 shows a comparison of the internal energy 
level in the aluminum shell material for the uncoated cylinder 
and the two composite cylinders, the latter exhibiting 
significantly higher values. The effects of the rubber 
coatings on the cylinder were significant and it is speculated 
that the adverse response of the coated cylinders is related 
to the energy build-up in the metallic material caused by the 
rubber coatings. 

In the next phase of study, the analysis focused on 
altering certain parameters to determine if they had any 
effect on cylinder response. In the second case, Composite 
Model 1 was used and the shear modulus of the rubber was 
reduced from 95.8 psi to 72.2 psi, with the coefficient of 
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COMPOSITE MODEL 1 


Figure 34: Deformation of (a) uncoated aluminum cylinder and 
(b) composite cylinder #1 at 2.66 msec.. 
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ALUMINUM MODEL 



Figur* 35: Defornacion of (a) uncoatec? aluminum cylinder and 
(b) composite cylinder 12 at 5.86 msec.. 
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Figure 36: Internal energy of aluminum shell material for 
uncoated and coated aluminum cylinders. 
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asymmetry held constant. Values for axial and hoop strain 
matched closely up to about 1.4 msec.. Thereafter, the 
magnitudes varied and the peaks occurred at different times. 
Plots of axial and hoop strain are shown in Figures 37-54. 
Again, the location was important in determining whether the 
effect of lowering shear modulus resulted in a larger 
response. Although the variation was small, lowering the shear 
modulus appeared to adversely affect the response. The 
internal energy of the aluminum material for the two dif'^erent 
values of shear modulus is shown in Figure 55. After 1 msec., 
the energy levels differed considerably. The cylinder with the 
lower shear modulus showed the higher internal energy. As in 
the previous study, the inertial effects on cylinde: response 
are most pronounced near the endplate ("C" element locations) . 

The shear modulus of the rubber was then increased with 
the coefficient of asymmetry held constant. As the modulus was 
increased by orders of magnitude, the response of the cylinder 
dramatically improved. Effective plastic strain is plotted for 
element A3 and is shown in Figure 56. For shear modulus values 
at or above 500 psi, the aluminum exhibited significantly 
lower strain. To help visualize how the response was affected 
by a change in the material property of the rubber, radial 
velocity at node 49 {interior node at element A3 along the 
symmetric boundary) was examined. The plot is shovm in Figure 
57 and it indicates that at values above 500 psi, the response 
is not erratic as in the low shear modulus case, but rather 
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Figure 37; Axial strain at position A1 for coated aluminum 
cylinders with variation of rubber shear modulus 
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Figure 30; Hoop strain at position A1 for coated aluminum 
cylinders, with variation of rubber shear modulus 
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Figure 39: Axial strain at position A2 for coated aluminum 
cylinders with variation of rubber shear modulus 
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Figure 41: Axial strain at position A3 for coated aluminum 
cylinders with variation of rubber shear modulus. 
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Figure 42: Hoop strain at position A3 for coated aluminum 

cylinders with variation of rubber shear modulus. 
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Figure 43: Axial strain at position B1 for coated aluminum 
cylinders with variation of rubber shear modulus 
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Figure 44: Hoop strain at position B1 for coated aluminum 
cylinders with variation of rubber shear modulus 































MICROSTRAIN 


AXIAL STRAIN FOR ELEMENT B2 



Figure 45: Axial strain at position B2 for coated aluminum 
cylinders with variation of rubber shear modulus. 
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Figure 46: Hoop strain at position B2 for coated aluminum 

cylinders with variation of rubber shear modulus. 
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'igure 47; Axial strain at position B3 for coated alurr.inuir. 

cylinders with variation of rubber shear modulus. 
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Flgurt 48: Hoop strain at position B3 for coated aluminum 

cylinders with variation of rubber shear modulus 
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Figure 49: Axial strain at position Cl for coated aluminum 
cylinders with variation of rubber shear modulus. 
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Figure 50: Hoop strain at position Cl for coated aluminum 

cylinders with variation of rubber shear modulus. 
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Figure 51: Axial strain at position C2 for coated aluminum 
cylinders with variation of rubber shear modulus, 
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Figure 52: Hoop strain at position C2 for coated aluminum 

cylinders with variation of rubber shear modulus. 
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Figure 53: Axial strain at position C3 for coated alurainum 
cylinders with variation of rubber shear modulus. 
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Figure 54; Hoop strain at position C3 for coated aluminum 

cylinders with variation of rubber shear modulus. 
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Figure 55: Internal energy of aluminum shell material of 

coated aluminum cylinders with variation of rubber 
shear modulus. 
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Figure 56: Effective plastic strain at position A3 for 

uncoated and coated aluminum cylinders with large 
variation of rubber shear modulus. 
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Figure 57: Radial velocity at node 49 for uncoated and coated 
aluminum cylinders with large variation of rubber 
shear modulus. 
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shows a gradual and smooth return to the steady state. The 
deformation of Composite Model 1 for different values of 
rubber shear modulus is shown in Figure 58. Larger 
deformation was visible in the cylinder that had a shear 
modulus of 95.8 psi.. However, for larger values of shear 
modulus, the amount of deformation is small and is 
commensurate with deformation observed in the uncoated 
aluminum cylinder model. The internal energy of the aluminum 
shell material is plotted for different values of rubber shear 
modulus and is shown in Figures 59 and 60. These plots clearly 
show that when the shear modulus is high, the response of the 
coated cylinder is almost identical tc the uncoated cylinder. 
The results of this study seem to indicate that the dynamic 
response of the cylinder is related to the stiffness of the 
coating and is therefore affected by the value of the shear 
modulus. Furthermore, the dramatic change in the response of 
the cylinder due to a change in the shear miOduius of the 
coating suggests the existence of a shear .modulus “threshold" 
value, i.e. a value below vchich causes t.he coated cylinder 
adversely respond to an underwater shock wave. 

In the third case, the ob 3 ective was to examine the 
sensitivity of the response to a change in the coefficient of 
asymmetry while holding the shear modulus of the rubber 
constant. Again, Composite Model 1 was utilized and "he 
coefficient of asymmetry was varied from 0.223 to 0.“50. At 
nearly all locations, axial and hoop strain values exhibited 
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COMPOSITE MODEL 1 (G=95.8 PSD 




ALUMINUM MODEL 



Figure 58; Deformation of uncoated and coated aluminum 

cylinders with large variation of rubber shear 
modulus at 2.86 msec. (Scale factor of 5). 
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Figure 59: Incernal energy of aluminum material for uncoated 
and coated cylinders with large variation of 
rubber shear modulus (10 msec, range). 
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Figure 60: Internal energy of aluminum material for uncoated 
and coated cylinders with large variation of rubber 
shear modulus (2 msec, range) 
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relatively small deviations until late time (later than 2 
msec.). The internal energy level of the aluminum, shown in 
Figure 61, was virtually identical for both cylinders up 
through 2.2 msec, and then showed a small deviation 
thereafter. Strain results are shown in Figures 62-79, and 
indicate that the dynamic response of the coated cylinder was 
wealcly influenced by changes in the coefficient of asymmetry 
and were not evident until late time. 


INTERNAL ENERGY OF ALUMINUM SHELL MATERIAL 
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Figure 61: Comparison of internal energy of aluminum shell 
material of coated cylinder with variation in 
coefficient of asymmetry. 
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Figure 62: Axial strain at position Al for coated aluminum 
cylinders with variation of coefficient of 
asymmetry. 
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Figure 63; Hoop strain at position Al for coated aluminum 
cylinders with variation of coefficient of 
asymmet ty . 
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Figure 64: Axial strain at position A2 for coated aluminum 
cylinders with variation of coefficient of 
asymmetry. 
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Figure 65: Hoop strain at position A2 for coated aluminum 
cylinders with variation of coefficient of 
asymmetry. 
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Figure 66: Axial strain at position A2 for coated aluminum 
cylinders with variation of coefficient of 
asymmetry. 
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Figure 67; Hoop strain at position A3 for coi'• ed v.luminum 
cylindeis with variation of costi.'.c;ent of 
asymmetry. 
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AXIAL STRAIN FOR ELEMENT B1 
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Figure 68: Axial strain at position B1 for coated aluminum 
cylinders with variation of coefficient of 
asymmetry. 
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Figure 69; Hoop strain at position B1 for coated aluminum 
cylinders with variation of coefficient of 
asymmetry. 
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Figure 70: Axial strain at position B2 for coated aluitinuin 
cylinders with variation of coefficient of 
asymmetry. 
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Figure 71: Hoop strain at position B2 for coated aluminum 
cylinders with variation of coefficient of 
asymmetry. 
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Figure 72: Axial strain at position B3 for coated aluxninuiT'. 
cylinders with variation of coefficient of 
asymr^etry. 
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Figure 73: Hoop strain at position B3 for coated aluminum 
cylinders with variation of coefficient of 
asjt'Ttmetry. 













































MICROSTRAIN 


AXIAL STRAIN FOR ELEMENT C1 
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Figure 74: Axial strain at position Cl for coated aluminum 
cylinders with variation of coefficient of 
asymmetry. 
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Figure 75: Hoop strain at position Cl for coated aluminum 
cylinders with variation of coefficient of 
asymmetiy. 
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Figure 76: Axial strain at position C2 for coated aluminum 
cylinders with variation of coefficient of 
asymmetry. 
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Figure 77: Hoop strain at position C2 for coated aluminum 
cylinders with variation of coefficient of 
asymmetry. 
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AXIAL STRAIN FOR ELEMENT C3 
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Figure 78: Axial strain at position C3 for coated aluminum 
cylinders with variation of coefficient of 
asymmetry. 
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Figure 79: Hoop strain at position C3 for coated aluminum 
cylinders with variation of coefficient of 
asymmetry'. 









In the fourth case, the objective was to determine how a 
variation in the thickness of the aluminum would affect the 


dynamic response of the coated cylinder. The thickness of the 
aluminum shell material in Composite Model 1 was changed from 
0.250 in. to 0.125 in. and 0.375 in.. The thickness of the 
aluminum endplate remained constant (1 in. thick). Plots of 
axial and hoop strain are shown in Figures 80-97. The results 
show that the thickness variation affected the cylinder 
differently depending on the location where an element was 
analyzed. For example, away from the endplate {"A" and "B" 
elements), hoop strain values were generally higher when the 
aluminum thickness was reduced, however the thinner aluminum 
resulted in lower axial strain values near the endlate. To 
gain a better perspective on what the overall effect was, 
effective plastic strain was plotted for ail elements cf 
interest and is shown in Figures 98-IC6. These plots showed 
that the thicker aluminum generally resulted in lower plastic 
strain at locations away from the endplate but higher plastic 
strain near the endplate. Alu-minum deformation is showr. in 
Figure 107 for the three different thickness cases. 
Deformations were scaled by a factor of five. It can not ho 
deduced which pattern of deformation represents t.he mast 
severe affect. The inter.nal energy of the aluminum shell 
material was also plotted and is shown in Figure 108. The 
results indicate that as the metal thickness was reduced, the 
internal energy was also reduced. Results cf this study seem 
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Figure 80: Axial strain at position A1 for coated cylinders 
with variation of aluminum shell thickness. 
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Figure 81: Koop strain at position Al for coated cylinders 
with variation of aluminum shell thickness. 
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AXIAL STRAIN FOR ELEMENT A2 
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Figure 82: Axial strain at position A2 for coated cylinders 
with variation of aluminum shell thickness. 
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Figure 83: 


TlME{sec) 

Hoop strain at position A2 for coated cylinders 
with variation of aluminum shell thickness. 
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Figure 84; Axiai strain at position A3 for coated cylinders 
with variation of aluitinuin shell thickness. 
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Figure 85; Hoop strain at position A3 for coated cylinders 
with variation of aluininum shell thickness. 
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Figure 86: Axial strain at position Bl for coated cylinders 
with variation of aluminum shell thickness. 


HOOP STRAIN FOR ELEMENT B1 


•0.9x10* 


— 0 375* THICK ALUMINUM 

- 0 250* THICK ALUMINUM 

•••• 0 125* THICK ALUMINUM 


TlME(sec) 

Figure 87: Hoop strain at position Bl for coated cylinders 
with variation of aluminum shell thickness. 
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?igxire 88s Axial strain at position B2 for coated cylinders 
with variation of aluminum shell thickness. 
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Figure 89s Hoop strain at position B2 for coated cylinders 
with variation of aluminum shell thickness. 
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Pigrure 90: Axial strain at position B3 for coated cylinders 
with variation of aluminum shell thickness. 
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Piflfure 91: Hoop strain at position B3 for coated cylinders 
with variation of aluminum shell thickness. 
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Figure 92: Axial strain at position Cl for coated cylinders 
with variation of aluminum shell thickness. 
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Figure 93: Hoop strain at position Cl for coated cylinders 
with variation of aluminium shell thickness. 
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Figure 94: Axial strain at position C2 for coated cylinders 
with variation of aluminum shell thickness. 
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Figure 95: Hoop strain at position C2 for coated cylinders 
with variation of aluminum shell thickness. 
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Figure 96; Axial strain at position C3 for coated cylinders 
with variation of aluminum shell thickness. 
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Figure 97: Hoop strain at position C3 for coated cylinders 
with variation of aluminum shell thickness. 
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EFFECTIVE PLASTIC STRAIN FOR ELEMENT A1 
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Figure 98: Effective plastic strain at position A1 for coated 
cylinders with variation of aluminum shell 
thickness. 
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Figure 99; Effective plastic strain at position A2 for coated 
cylinders with variation of aluminum shell 
thickness. 
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EFFECTIVE PLASTIC STRAIN FOR ELEMENT A3 
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EFFECTIVE PLASTIC STRAIN FOR ELEMENT B2 
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Figure 102: Effective plastic strain at position B2 for coated 
cylinders with variation of aluminum shell 
thickness. 
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Figure 103; Effective plastic strain at position B3 for coated 
cylinders with variation of aluminum shell 
thickness. 
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Figure 104; Effective plastic strain at position Cl for coated 
cylinders with variation of aluminum shell 
thickness. 
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Figure 105: Effective plastic strain at position C2 for coated 
cylinders with variation of aluminum shell 
thickness. 























MICROSTRAIN 


EFFECTIVE PLASTIC STRAIN FOR ELEMENT C3 



Figure 106: L;feetive plast 
cylir.iei's -.v^th 
chi;:kne3s. 


TIME(S^) 

ic stra:.n at. position C3 for coated 
vdMd11 on of dAurr\in uiTi s h6i 1 


78 





















0.375* THICK ALUXINoM (DISPL. SCALE FACTCK OF 5> 

Flgurfi 107: Defonnaticn of coated cylinders with variation of 
aluff.inuin shell thickness at 2.34 msec.. 
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Figuxt 108j Internal energy of aluininum shell material of 

coated cylinders with variation of aluminum shell 
thickness. 
























































to indicate that variations of metal thickness have a local 
effect on response and are greatly influenced by inertia 
effects. 

In the fifth case, the objective was to determine how a 
variation in the coating thickness would affect the dynamic 
response of the cylinder. The thickness of the tread stock 
rubber in Composite Model 1 was changed from 0.250 in. to 
0.125 in. and to 0.375 in.. Axial and hoop strain values are 
plotted and are shown in Figures 109-126. At nearly all 
locations, the response of the cylinder significantly improved 
when the rubber thickness was increased but was adversely 
affected when the rubber thickness was decreased. Early 
termination of the proble.m in the thin rubber case was due to 
a severe distortion of the model mesh. The effects of rubber 
thickness are also evident in the values of effective plastic 
strain that were achieved and are shown in Figures 127-135. 
Differences in response can be seen in the deformation of the 
cylinder for each case and are shown in Figures 136 and 137. 
The effect is more dramatic at lace time and can be seen more 
clearly in the response at 2 86 msec. (Figure 137). Cylinder 
response was also examined for a rubber thickness of 0.500 in. 
and it was found that the response of the cylinder improved. 
Inte.rnal energy of the aluminum shell material is plotted for 
the various rubber thicknesses and is shown in Figure 138. As 
observed in case 2 (shear modulus variation), the results of 
this study also support the hypothesis that the dynamic 
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Figure 111: Axial strain at position A2 for coated aluininum 
cylinders with variation of rubber thickness. 
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Figure 112: Hoop strain at position A2 for coated aluminum 
cylinders with variation of rubber thickness. 
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AXIAL STRAIN FOR ELEMENT A3 
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Figure 113: Axial strain at position A3 for coated aluminum 
cylinders with variation of rubber thickness. 
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Figure 114: Hoop strain at position A3 for coated aluminum 
cylinders with variation of rubber thickness. 
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MICROSTRAIN 


AXIAL STRAIN FOR ELEMENT B1 



TlME(sec) 

Figure 115; Axial strain at position B1 for coated aluminum 
cylinders with variation of rubber thickness. 
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Figure 116: Hoop strain at position B1 for coated aluminum 
cylinders with variation rubber thickness. 






















































MCI 


Figure 117: Axial strain at position B2 for coated aluminuiT^ 
cylinders with variation of rubber thickness. 
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Figure 118: Hoop strain at position 82 for coated aluminum 
cylinders with variation of rubber thickness. 
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- 0 375* THICK RUBBER 

- 0 250* THICK RUBBER 

-0.125* THICK RUBBER 


TlME(sec) 

figure 119: Axial strain at position B3 for coated aluminum 
cylinders with variation of rubber thickness. 
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Figure 120: Hoop strain at position B3 for coated aluminum 
cylinders with variation of rubber thickness. 
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Pi(jurtt 121: Axial strain at position Cl for coated aluminum 
cylinders with variation of rubber thickness. 
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Figure 122: Hoop strain at position Cl for coated aluminum 
cylinders with variation of rubber thickness. 
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Figure 123: Axial strain at position C2 for coated aluminum 
cylinders with variation of rubber thickness. 
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Figure 124: Hoop strain at position C2 for coated aluminum 
cylinders with variation of rubber thickness. 
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AXIAL STRAIN FOR ELEMENT C3 
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Figure 125: Axial strain at position C3 for coated aluminum 
cylinders with variation of rubber thickness. 
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Figure 126: Hoop strain at t^s^ition C3 for coated aluminum 
cylinders with variation of rubber thickness. 
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EFFECTIVE FUSTIC STRAIN FOR ELEMENT A1 


0.375- THICK RUSTIER 
0.250* THICK RUBBER 
0.125* THICK RUBBER 


TIMEtSW) 

Flgfure 127: Effective plastic strain at position A1 for coaved 
aluminum cylinders with variation of rubber 
thickness. 
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Figur* 126: Effective plastic strain at position A2 for coated 
aluminum cylinders with variation of rubber 
thickness. 
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Figure 129: Effective plastic strain at position A3 for coated 
aluminura cylinders with variation of rubber 
thickness. 
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EFFECTIVE PLASTIC STRAIN FOR ELEMENT B2 
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Pijjure 131: Effective plastic strain at position B2 for coated 
alumin'om cylinders with variation of rubber 
thickness. 
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Figurt 132: Effective plastic strain at position B3 for coated 
aluminuit cylinders with variation of rubber 
thickness. 

































MJCROSTRAIN 


EFFECTIVE PLASTIC STRAIN FOR ELEMENT Cl 



- 0.375-THICK RUBBER 

- 0.250* THICK RUBBER 

-0.125* THICK RUBBER 


TIME(sec) 

Figure 133: Effective plastic strain at position Cl for coated 
aluir.inum cylinders with variation of rubber 
thickness. 
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Figure 134: Effective plastic strain at position C2 for coated 
aluminum cylinders with variation of rubber 
thickness. 








































































MICROSTRAIN 


EFFECTIVE PLASTIC STRAIN FOR ELEMENT C3 
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Figure 135; Effective plastic strain at position C3 for coated 
aluminum cylinders with variation of rubber 
thickness. 
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Figur# 136: Deformation of coated aluminum cylinders with 
variation of rubber thickness at 0.76 msec. 
(Displacement scale factor of 5). 
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Flgurs 137: Deformation of coated aluminum cylinders with 
variation of rubber thickness at 2.86 msec. 
(Displacement scale factor of 5). 
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response of the cylinder is related to the stiffness of the 
coating and is therefore affected by the thickness of the 
coating. The significant difference in response due to changes 
in coating thickness again suggest the existence of a 
thickness "threshold" value, i.e. a value below which causes 
the coated cylinder to respond adversely to an underwater 
shock wave. 

In the sixth case, the objective was to determine how a 
change in the metal material would affect the dynamic response 
of the cylinder. For this study, an uncoated steel cylinder 
was compared to a steel cylinder coated with tread stock 
rubber. The steel shell thickness and coating thickness were 
both 0.250 in.. The endplate thickness was 1 in. as in the 
aluminum cylinder cases. Axial and hoop strain values are 
plotted for all element locations and are shown in Figures 
139-156. Results appeared analagous to those obtained in case 
1. Again, axial and hoop strain values were larger in the 
coated cylinder than in the uncoated cylinder. It was noted 
that the ratio of the maximum response of the coated and 
uncoated cylinders was generally higher in the case of steel 
than in the case of aluminum. For example at location B3, for 
response up through 1.3 msec., the ratio of coated to uncoated 
maximum axial strain was 4.3 for steel and 1.8 for alum.inum. 
This is likely attributable to the strain-rate sensitivity of 
steel. Results of this study indicate that effects of surface 
coatings on the response of cylinders subjected to underwater 
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Figure 139: Axial strain at position A1 for uncoated and 
coated steel cylinders. 
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Figure 140: Hoop strain at position A1 for uncoated and coated 
steel cylinders. 
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figure 141; Axial strain at position A2 for uncoated and 
coated steel cylinders. 
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Figure 142: Hoop strain at position A2 for uncoated and coated 
steel cylinders. 
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Axial strain at position A3 for uncoated and 
coated steel cylinders. 
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Figure 144: Hoop strain at position A3 for uncoated and coated 
steel cylinders. 
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Figure 145: Axial strain at position B1 for uncoated and 
coated steel cylinders. 
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Figure 146: Hoop strain at position B1 for uncoated and coated 
steel cylinders. 
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Figure 147; Axial strain at position 
coated steel cylinders. 


B2 for uncoated and 
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Figure 148: Hoop strain at position B2 for uncoated and coated 
steel cylinders. 
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Figure 149: Axial strain at position B3 for uncoaced and 
coated steel cylinders. 
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Figure 150: Hoop strain at position B3 for uncoated and coated 
steel cylinders. 








































rflCROSTRAIN _ MICROS i RAIN 


AXIAL STRAIN FOR ELEMENT C1 



Figure 151i 


TIME(sec) 

Axial strain at position Cl for uncoated and 
coated steel cylinders. 
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Figure 152; Hoop strain at position Cl for uncoated and coated 
steel cylinders. 


106 




































































MICROSTRAIN 


AXIAL STRAIN FOR ELEMENT C2 



[11 


























































MICnOSTRAIN 



























































shock are similar, and while the magnitudes of the response 
vary, the trends are the same. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


Under certain conditions, when subjected to underwater 
shock, compliant coatings appeared to concentrate wave energy 
within the structure for longer duration, resulting in 
significantly higher magnitudes of stress and strain. Instead 
of a gradual release of energy into the surrounding water 
medium, most of the energy was retained in the metal material. 
Cylinder response was most influenced by changes in rubber 
shear modulus and rubber thickness. Increasing these values 
resulted in improved cylinder response. Both parameters are 
related to the coating stiffness and it is this property that 
most likely governs the extent of the energy transfer in the 
compliant coating and consequently the dynamic response of the 
cylinder. Results point to the existence of threshold values 
for both coating shear modulus and coating thickness, values 
below which lead to adverse and erratic cylinder response. 

A more detailed analysis of wave transmission and 
reflection at the rubber-metal and rubber-water interfaces 
would be beneficial. An examination of the effects of 
stiffness changes at the microscopic level would improve our 
understanding of the energy transfer mechanism in the coating. 
A sensitivity analysis should also be conducted to investigate 
cylinder response for changes in endplate thickness and 
loading (asymmetrical vice symmetrical). 


no 








Follow on physical tests should be conducted to validate 
the numerical models and tc- enhance our understanding of the 
effects of compliant coatings on shock hardening. 
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APPENDIX - COMPOSITE CYLINDER MODEL INQRID INPUT FILE 


MODEL-MOONRIV,KING //THICKSHELL AND BRICK, ALUMINUM AND RUaBER(TREAD STOCK)// 
dn3d vec term 0.010 plti 20.0e-6 prti 1000.0 


mat 1 type 3 e 1.08e+7 pr 0.33 ro 2.610e-4 etan 0.0 sigy 4.0o+4 
tshell quad 3 thick 0.250 endmat 

mat 2 type 3 e 1.08e’(>7 pr 0.33 ro 2.160e-4 etan 0.0 sigy 4.0e'('4 
tshell quad 3 thick 1.000 endmat 

mat 3 type 3 e 1.08e+7 pr 0.33 ro 2.160e-4 etan 0.0 sigy 4.0e+1 
endmat 

mat 4 type 3 e 1.08e+7 pr 0.33 ro 2.160e-5 etan 0.0 sigy 4.0e+4 
endmat 


mat 5 type 3 e 1.08e+7 pr 0.33 ro 2.160e-5 etan 0.0 sigy 4.0e+4 
endmat 


c Materials 3''5 are modeled as the rubber material, however the Ingrid program 
c will not accept this material specification. Therefore, these materials will 
c be entered as if aluminum brick elements and the correct properties for the 
c rubber will be edited in the ingrido file. 


c The material input line for the rubber material parts should appear as 
c follows: 

c mat 3 type 27 a 29.291 b 18.609 pr 0.4991 ro 9«20e-05 endmat 

c mat 4 type 27 a 29.291 b 18.609 pr 0.4991 ro 9.20e-05 endmat 

c mat 5 type 27 a 29.291 b 18.609 pr 0.4991 ro 9.20e-05 endmat 

c Materials 1 and 2 are Aluminum and are used for the inner cylindrical shell 

c and inner end plate. Materials 3,4 and 5 are Tread Stock (Type of rubber) 
c and are used for the outer cylindrical shell and outer end plate. 


led 1 2 0.0 

0,0 

5.0e-2 

0.0 




led 2 2 0.0 

2.0e-7 

5.0e-2 

2.0e-7 




plan 2 

0 0 0 0 1 0 

0.001 

symm 





0 0 0 1 0 0 

0.001 

symm 





c name 

type 

point 

dir vector 

size 



sd 1 

cyll 

0 0 0 

0 1 0 

6.250 



sd 2 

cyli 

0 0 0 

0 1 0 

6.000 



8d 3 

cyli 

0 0 0 

0 1 0 

5.750 



sd 4 

cyll 

0 0 0 

0 1 0 (5 

. 750*2./3.] 



c ********** 

Generate Inner Cylindrical Shall ******.... 


start 







1 2 

3 

6 

9 10 

U > 



1 2 

4 

6 

6 13 

18 20 22 

24 

23 ) 

1 2 

3 

6 

9 10 

11 ? 



-1. -1 

. -1. 

0. 

1. Is 

1. 



-21.0 -20 

.0 -17. 

5 -14.5 -11.0 0.0 

11.0 14.5 17.5 

20.0 

21.0 

1. 1 

. 1. 

0. 

-i. -i. 

-la 

‘ 



or 600700kl 
or lOOaOOJk 
or 0 0 1 0 0 2 i J 
or 006007Ji 
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di 13057; J13057; 
d 2 0 2 6 0 6 

sfi -1 -7 ; ; -1 -7 ; sd 2 
sfl -2 -6 ; ; -2 -6 ; sd 3 

« dOlOOGO 

d 1 0 0 4 0 0 

mate 1 

* 

end 


c ********** Generate Inner End Plate ********** 
start 


1 

3 6 

9 

11 ; 

1 

2 ; 



1 

3 6 

9 

11 ; 

-1. 

-1. 0. 

1. 

1. 

20.0 

21.0 



-1. 

-1. 0. 

1. 

1. 

or 0 0 0 

0 0 0 k i 



di 1 2 

0 4 5 ; ; 1 2 

0 4 5 

• 

/ 

sfl -1 

-5 ; ; -1 -5 

; sd 

3 

sfvi -2 

-4 ; ; -2 -4 

; sd 

4 

d 1 0 0 

3 0 0 




mate 2 


end 

c ********** Generate Outer Cylindrical Shell Coating ********** 
start 

1 2 3 6 9 10 11 ; 

\ I 3 5 7 9 14 19 21 23 25 26 27 ; 

1 2 3 6 9 10 11 ; 

“I* “!• “1. 0. 1. 1, 1, 

-21.25 -21.0 -20.0 -17.5 -14.5 -11.0 0.0 11.0 14.5 17.5 20.0 21.0 21.25 

!• 1 . 0 . - 1 . - 1 . - 1 . 

or 600700kj 
or 1 0 0 2 0 0 J k 
or 001002ij 
or 006007ji 

pri -1 -7 j ;-l -7 ; 1 -1.0 0.0 0.0 0.0 

dil3057;;13057; 
d 2 0 2 6 0 6 

sfi -1 -7 ; j -1 -7 ; sd 1 
sfi -2 -6 } ; -2 -6 ; sd 2 

d 0 1 0 0 7 0 
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d 1 0 0 4 0 0 


mate 3 
end 

c ********** Generate Outer End Plate Coating (Small Ring Area) 


start 



1 

2 

3 

6 

9 

10 

11 ; 


1 

2 ; 







1 

2 

3 

6 

9 

10 

11 ; 


-1. 

-1. 

-1. 

0. 

1. 

1. 

1. 


21.0 

21.250 







-1. 

-1. 

-1. 

0. 

1. 

1. 

1. 

or 

6 0 0 

7 0 0 k 

j 





or 

1 0 0 

2 0 0 j 

k 





or 

0 0 1 

0 0 2 i 

j 





or 

0 0 6 

0 0 7 j 

i 






dil3057;;13057; 
d 2 0 2 6 0 6 


sfi -1 -7 ; ; -1 -7 ; sd 2 
sfi -2 -6 ? ; -2 -6 ; sd 3 

d 1 0 0 4 0 0 


mate 4 


end 

c ********** Generate Outer End Plate Coating (Large Plate Area) 


start 


1 

3 6 9 

11 ; 

1 

2 ; 


1 

3 6 9 

11 ; 

-1. 

-1. 0. 1. 

1. 

21.0 

21.250 


-1. 

-1. 0. 1. 

1. 

or 0 0 0 

0 0 0 k i 


pr 1 2 1 

525 1 -1.0 

0.0 0 

di 1 2 

0 4 5 ; ; 1 2 0 

4 5 

sfi -1 

-5 ) ) -1 -5 ; 

sd 3 

sfvi -2 

-4 J > -2 -4 ; 

sd 4 

d 1 0 0 

3 0 0 


mate 5 




end 


end 


********** 


**★*****★* 
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